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We consider a configuration consisting of a wormliole filled by a perfect fluid. Such a model can 
be applied to describe stars as well as neutron stars with a nontrivial topology. The presence of a 
tunnel allows for motion of the fluid, including oscillations near the core of the system. Choosing 
the polytropic equation of state for the perfect fluid, we obtain static regular solutions. Based on 
these solutions, we consider small radial oscillations of the configuration and show that the solutions 
L~:> ' are stable with respect to linear perturbations in the external region. 
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I. INTRODUCTION 



After the discovery of the accelerated expansion of the present Universe at the end of the 1990ies, a great deal of 
attention has focused on models in which violation of one of the energy conditions takes place. In hydrodynamical 
language, this means that the parameter if in the equation of state p = wp, defining the ratio of the pressure p of 
some matter to its energy density, may be less than —1/3 (i.e. violation of the strong energy condition), or even less 
^ • than —1 (i.e. violation of the weak energy or null energy condition). While w = — 1 is realized by the presence of 
I , a cosmological constant, leading to exponential expansion of the Universe, w < —1 would be realized by phantom 
2 • matter with a possible end of the Universe in a finite time. Recent astronomical observations indicate that such a 
> \ possibility is not excluded P, Q • 
. On the other hand, when considering models of compact astrophysical objects, a violation of the weak/null energy 
'—i' condition leads to a possible existence of configurations with non-trivial topology - traversable wormholes Q. At 
the present time, much work has been devoted to the study of various models of both microscopic and macroscopic 
' wormholes. One of the possible variants of obtaining wormhole-like solutions consists in considering ghost scalar fields, 
', i.e. scalar fields with the opposite sign in front of the kinetic energy term. The use of such fields allows, in some cases, 
■ to obtain an equation of state with u; < — 1, and thus to violate the weak/null energy condition. 
\ Models with ghost scalar fields have been considered earlier. The early pioneering works with a massless ghost scalar 
■"slj" . field have been done in Refs. [1, @. A ghost scalar field with a Mexican hat potential was investigated in the papers 
' @i01i where it was found that this system had regular, stable solutions only for topologically non-trivial (wormhole- 
. like) geometry. In (8l-[l3| traversable Lorentzian wormholes were investigated further, refining the conditions on the 
type of matter or fields that would lead to such space-times. A general overview of the subject of Lorentzian wormholes 
and violations of the various energy conditions can be found in the book by Visser [T5| . 

In the case of macroscopic wormholes, ghost fields can be used to create models of astrophysical objects which can 
serve as entrances to wormholes. From the point of view of a distant observer, such objects will be quite similar 
to the usual star-like configurations, but they will have some characteristic features. Moreover, in the presence of 
electric and/or magnetic charges there exist new possibilities to detect such objects via galactic and cxtragalactic 
observations. A discussion on the observability of such electrically /magnetically charged wormholes is given in [l6j . 

However, it is also possible to imagine a situation where the concentration of exotic matter violating the weak/null 
energy condition is formed at the center of an otherwise ordinary star. This in turn then leads to the possibility of 
the creation of a tunnel - a wormhole whose throat is filled with ordinary (star) matter. Thus an ordinary star or a 
neutron star could harbour some exotic matter at its core providing a nontrivial wormhole topology. Such a combined 
configuration with a nontrivial spacetime topology differs from the usual spherically symmetric stars consisting of 
ordinary (nonexotic) matter and having a center at the point r = 0. In the model considered here, there is no such 
center but there exists some minimal value of the radial coordinate r — r„iin corresponding to the radius of the 
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throat. But such a configuration will still be different from a usual wormhole due to the presence of ordinary matter. 
Taking this into account, we will describe the region around the throat by using the term "core" along with the term 
"throat" . The advantage of introducing this term is that on the one hand it implies the presence of the throat and, on 
the other hand, it allows to keep the spirit of considering star-like configurations consisting only of ordinary matter. 

For a distant observer, such a configuration would very much look like an ordinary star. However, also some 
distinctions should be present. For example: 

(i) One of the most striking differences is that instead of a single star two similar stars will be observed that are 
separated in space. These two stars are associated with the two mouths of the wormhole. 

(ii) Due to the fact that matter inside such star can move freely through the tunnel, the presence of radial (including 
quasiperiodical) motion is possible. Thus one can estimate the energy associated with such motion. 

In this connection, we here suggest to consider a simple (toy) model for a "wormhole plus ordinary matter" , where 
the source of the exotic matter is a massless scalar field, and the ordinary matter is a perfect fluid with a polytropic 
equation of state p = kp^ . In section HIl we consider the general properties of this model. In section Hill we find 
static solutions for the perfect polytropic fluid filling the wormhole's throat. In section IIVI the infinitesimal radial 
oscillations of the system are considered. In Appendix |A] the equations for the perfect polytropic fluid are derived 
following standard notation. Finally, in Appendix [B] the oscillating solutions for the perfect fluid with a polytropic 
equation of state are found; while in Appendix C stability of the external solution under linear perturbations is 
addressed. 



II. GENERAL PROPERTIES OF THE MODEL 



Let us first consider some general properties of the model consisting of a gravitating massless ghost scalar field 4> 
with negative kinetic term and a normal isotropic perfect fluid. The Lagrangian of the system is chosen in the form 

L = -^-ia,0a^<^ + L™, (1) 

where Lm refers to the isotropic perfect fluid. The energy-momentum tensor of this system is then 

T^ = {p + p)u,u'^ - 5^p - a,^^'^^ + ^<5f a^</)9''0, (2) 

where p and p are the energy density and the pressure of the fluid, u' is the four- velocity (in units where c = 1). 
For the metric of this system we employ Schwarzschild-like coordinates 

ds^ = e'^^^'Ue - (de^ + sin^ . (3) 

1 — b[r)/r ^ ' 

The general properties of such a metric describing a wormhole can be found in Refs. (sl. fisl. [TH : 

(i) the coordinate r covers the range <.r < -t-oo, where rg is the throat radius; 

(ii) in order to cover the whole space-time one needs to use two copies of the coordinate system ([3]); 

(iii) the requirement of the absence of horizons or singularities leads to the condition that v should be everywhere 
finite; 

(iv) throughout space-time (1 — h(r)lr) > implying 

b{ro) = ro, 6'(ro) < 1, 6(r) < r; 

(v) to provide asymptotic flatness of the system, it is necessary to satisfy the asymptotic limit 

b{r)/r — > as |r| — > oo. 

Taking into account these general properties of the model under consideration, in studying the "wormhole plus 
star" system it is necessary, first of all, to choose some equation of state for the fluid filling the wormhole. As an 
example, we next consider a perfect fluid with a polytropic equation of state, which is of a considerable interest in 
astrophysical applications. 
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III. PERFECT FLUID WITH THE EQUATION OF STATE p = k p"' 



We now consider the case that the fluid filling the wormhole has a polytropic equation of state. In our calculations, 
we will follow Ref. [l7| . where the model of a general relativistic polytropic fluid sphere is considered. The equation 
of state for the polytropic fluid is 



p = kp^, 



(4) 



where fc, 7 are constants, and p is the mass density of the fluid. (The mass density corresponds to the energy density 
because of c = 1.) Follow Ref. (O], let us rewrite this equation in terms of the variable 9 introduced as 



P = PcO" 



(5) 



where pc is the core density, and the constant n, called the polytropic index, is related to 7 via n ~ 1/(7 ~ !)• Then 
Eq. (HJ takes the form 



Using this equation of state, one finds the following system of equations (for details, see Appendix \K\ : 
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d^ ^ 
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-2(n+iy 



(6) 

(7) 

(8) 



where D is the dimensionless constant 



^ is the dimensionless radius 



and i'(^) is the dimensionless function 



D 



AttGD 
(n + l)cr' 



A = 



A-kGpc 



-1 1/2 



v[Q = u r). 



(9) 



(10) 



The function uir) is a new metric function determining the component grr [see Appendix El Eq. (|A.8|) ]. and M is the 
total mass of the configuration within the range ^0 ^ ^ ^ Cb; where ^0 is the throat radius, and ^f, is the boundary of 
the fluid where = 0. The constant a is defined by the formula (jA.ll[) . and is equal to the ratio of the core pressure 
Pc and the core density pc of the fluid. 

Equations ([7]) and ^ are to be solved for given n and a subject to the boundary conditions 



9o = 0(?o) = l, fo = t'(eo) 



2(n + l)fT 



^0. 



(11) 



The first condition corresponds to the fact that the core density is pc, and the second condition follows from the fact 
that at the throat fo(ro) = rp (see condition (iv) of the previous section). Then, by using the definition of h{r) through 
u{r) (see Eq. (|X8|) 



h{r) = 2GMu{r) 



(12) 



one obtains condition (|TT|) . 

The boundary condition for w(^o) implies that the coefficient (1 — 2a {n + l)v/^) in front of the derivative of the 
function 9 in Eq. ([7]) goes to zero at the throat, which, in general, leads to the occurrence of a singularity in the 



solution. To avoid this, expand the solution of the system ©-(El) in the neighborhood of the point C = Co and start 
calculating the solution at ^ = choosing the boundary condition for v in the form 



Substituting this into Eq. ([5]), wc find the following expression for vi: 



Vl 



ID' 



1+a 

1+<T0O 



-2(ti+1) 



(13) 



(14) 



Further, to exclude the appearance of a singularity in Eq. ([7]), it is necessary to assume that, together with (1 — 2a{n + 
l)w/C), the right hand side of Eq. ([7]) simultaneously goes to zero. Proceeding from this requirement, we obtain the 
following expression for D: 



1 



2(n+l) 



1 



^0 



Substituting this into Eq. p4)) . we find 



(15) 



Thus, the mathematical description of an equilibrium static configuration is achieved by means of the system of 
equations ©-(jS]) with the boundary conditions Eq. (jlip for 6*0, and Eqs. ([T5|) and for w(Ci)- 

Examples of numerical solutions obtained in this way are presented in Figs. [T] and [5] for the polytropic index n ~ 1.5 
and n = 2, respectively. In these figures, the following functions are shown: the metric functions gtt = e'^ and 



Qrr = — e ; the distribution of the energy density of the fluid p, and the total energy density e = 
consisting of the scalar field and the fluid. (Both p and e are given in units of the core density pc.) 
For the graphs, we have employed the following expressions: e'^ is obtained from Eq. (jA.12[) 

2(ji+l) 



of the system 



(16) 



follows from Eq. (|A.8|) . using © and (dU]) 



(17) 

(note that this expression diverges at ^ = as discussed in connection with the boundary conditions), and e follows 
by using ((5|). ((9)) and (|T6)) 

The integration constant i>c, corresponding to the value of ly at the throat, is determined by requiring e'' to be equal 
to unity at infinity, corresponding to an asymptotically flat space-time. To fix the constant, it is necessary to find the 
external solution of the field equations outside the fiuid. For this purpose, let us use the Einstein equations ()A.5|) and 
(jA.6p . and the first integral (|A.4p . taking into account that in this region there is no ordinary matter, i.e. 9 — 0. This 
leads to the following system of equations: 



l-2cr(n + l) 



1 D'e-"" 



1 + 0- 



-2(ti+1) 



(18) 



5 




FIG. 1: The metric functions gtt = e"^, loglgrri ~ A, tiie fluid density p, and the total energy density e from psp for the 
polytropic index n = 1.5. The point where p = corresponding to the boundary of the fluid is situated at ^f, « 1.79. To 
provide the asymptotical flatness of the solutions, i.e. e'^,e^ — >■ 1 as ^ — ^ oo, the value of the constant Uc should be chosen as 
i^c ~ -1.46. 
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FIG. 2: The metric functions gtt = e", logjgrri ~ A, the fluid density p, and the total energy density e from p8|) for the 
polytropic index n = 2.0. The point where p = corresponding to the boundary of the fluid is situated at ^f, 1.99. To 
provide the asymptotical flatness of the solutions, i.e. e'^,e^ ^ 1 as ^ — > oo, the value of the constant Vc should be chosen as 
i^c ~ -1.57. 



which, by analogy with the transformations made above, can be rewritten in terms of the dimensionless variables 
v{^),(p{£,) and as follows 



dtp 



dv 
d^ 

dv 
d^ 

2 



'2^ 



1 



l-(7(n + l)|re"'" 
1 - 2cr(n + 1)t 



^4 1- 2(7(71 +1)|' 



(22) 
(23) 
(24) 



This system contains the parameter ct as a trace of the influence of the fluid on the external solution. The solution 
is to be sought beginning from the surface of the star at ^ = ^6 by using, as the boundary conditions, the values of 
v{£,b), 'fii^b) and vi^b) obtained from the solution of the equations ([7]) and © for the internal part of the configuration. 
This then allows to determine the value of the integration constant Vc by requiring e"^ to be equal to unity at infinity, 
providing asymptotical flatness of the space-time. (The values of for the examples shown in figures [1] and [5] are 
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given in the captions.) 

The systems of equations dZ])-® and ([^ - ([M]) can be used to determine the total mass of the system under 
consideration. For the spherically symmetric metric (jSj the function b{r) is associated with the effective mass m(r) 
inside the radius r as follows 0, lisl. Il5| 

m{r) = b{r)/2 = ro/2 + An f e{x)x^dx. 

J ra 

Here the integration constant is chosen to provide &(ro) = tq (see condition (iv) from section , and e(x) is taken 
from Eq. (|18l) . Using formulas ((TU)) and ((T^ . we introduce the dimcnsionlcss function 

B{^) = Ab = 2a{n+l)v{^) 

defining the dimensionless mass function of the configuration as M(^) = B{^)/2. Then the asymptotic value 
limj_>oo -^(C) = corresponds to the total mass of the configuration. For the cases presented in figures [1] and [H 
the total masses are Mn=i.^ ~ A/n=2.o ~ 0.37. Thus despite the violation of the energy conditions and the presence 
of a negative energy density, the total mass of the configurations remains positive. 



IV. RADIAL OSCILLATIONS OF THE SYSTEM 



We now study the radial oscillations of the above static system. Bearing in mind the possibility of using an arbitrary 
equation of state p = p{p)^ we will derive perturbation equations without imposing any preliminary restrictions on 
the equation of state. Various general-relativistic stellar models have been thoroughly investigated. In particular, 
infinitesimal radial oscillations of a gas sphere with polytropic equation of state has been discussed by Chandrasekhar 
[l8| . Here we employ the techniques developed in 18 1 for investigating the infinitesimal radial oscillations of the fluid 
filling the wormhole. 

To this end, it is more convenient to choose the metric in the form (following Ref. [HI): 



ds^ = e^dt^ - e^dr^ - 



sm 



(25) 



where the metric functions ly, X depend both on time t and radial coordinate r. The corresponding Einstein equations 
can be written in the form 
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■ 1 


A'' 


1 
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v'' 
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+ — 
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_}_d_ 



,-Al 



1 



(26) 
(27) 
(28) 



In the above equations, "prime" and "dot" denote differentiation with respect to r and t, respectively. Considering 
only infinitesimal oscillations the components of the four- velocity are given by 



with three- velocity 



dr 



and the index of the metric functions indicates the static background solutions of the Einstein equations. The 
components of the energy-momentum tensor ([5]) are then given by 



Tl 



-p+-e 



-A i,/2 



^0 = (po + Po)v + e ■ 



(29) 

(30) 
(31) 
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Let us seek perturbed solutions of the form 



y = yo + Vp, 



(32) 



where the index refers to static solutions, the index p indicates the perturbation, and y denotes one of the functions 
X,i^,p,p or (j). Substituting these expressions into Eqs. (p6|) and ((27)) . we find: 



-,-Ao 



rX'p + Ap (1 - rX'„)] = ^ [re-^"Xp] = 8^G 
A, 



r 



— Xr, 



8ttG 



-Pp 



'0 1 Vp 



:4>Q>^p 



Next, from Eq. (pSj) we have 



-Ao 



-Ji = -8ttG 
r 



The Einstein field equations arc not all independent since 
The i = I component of this equation has the form 



1/ .X 1 



In the model under consideration the perturbed components r| and Tg are given by 



= n 



-pp 



-Ac J,' A' 



\Vp 



^<t>'o>'p 



Now it is convenient to introduce a "Lagrangian displacement" C with respect to time t defined by 

Then Eq. ([55]) can directly be integrated to give 

Xp - -SttGc^V [{po +Po)C + e-^"^'o^p] . 
Using this expression, we have from Eq. (|33p 

Pp = e-^o^'o (^0; - ^4>o^p^ [(Po+Po) C + e-^"^'o^p] }, 



and the corresponding expression for i/'p follows from Eq. (p4l 

A, 



' r 



-Pp 



-Ao J,' 



Now we suppose that all perturbations have the following harmonic dependence on time t 

ypit,r)=yp{r)e^'^\ 



(33) 
(34) 

(35) 



(36) 



(37) 



(38) 



(39) 



(40) 



(41) 



where the function yp(r) depends only on the space coordinate r, and the characteristic frequency lo should be 
determined from the calculations. (For convenience, we hereafter drop the bar.) Then, using the expressions (|29p - (PT|) 
and (|37p . the perturbed equation ((36| takes the form 



^2gAo-..o 







(4>p - ^<P'o\^ - 


f0o 



{pp + Pp) v'q + e 



-An 



'I^Wo4>'p + \(t>'^ Wp - ^nK) + i^'p - \4>'o>^P^ - \{po+ Po) v'p = 0- 



(42) 
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By adding to this equation the equation for the perturbed scalar field <f>p which follows from 

1 d 



-0, 



and takes the form 



(43) 



we have a system of two second-order ordinary differential equations ([^^ and with respect to the scalar field 
perturbation (f>p and displacement C, which describes the infinitesimal radial oscillations of the system "wormhole plus 
star" . 

Thus, substituting equations ((38|) -(|40 )) . (|B.1[) and (jB.2p to the system (|42| and p3)) . regw/ar solutions of the indicated 
system are to be sought. As shown in previous investigations of such a ty pe o f problem, the solution of an equation 
similar to Eq. even without a scalar field, is a non-trivial problem [1^. In the latter case, when considering 

an ordinary spherically symmetric star-like configuration (where = 0), the boundary conditions are chosen in the 
following form [isj : 



C = at 



and = at r ~ R, 



where R is the radius of a star. As pointed out in [l9|, in the case of a polytropic equation of state which is widely 
used in the construction of hot (T > 0) stellar models of physical interest, 7 is finite and non-zero at r = i?, and po/po 
has a pole of order 1. In such a situation the point r = i? is a regular singular point of the eigenequation (|42l) . and 
one should choose the value of the derivative {d(/dr)r=B, to provide regularity of the solution at r = i?. The further 
solution then reduces to employing one of the methods described in Ref. [l^ . 

The numerical calculations for the case of a polytropic equation of state presented in Appendix |B] show that in the 
presence of a wormhole the situation on the boundary of the fluid remains similar to the case of an ordinary star: at 
r = R there is the divergence of 



CONCLUSIONS AND REMARKS 



The main purpose of the paper has been the demonstration that regular static solutions exist for the system of a 
star with a wormhole, describing self-consistently a configuration composed of a gravitating perfect fluid and exotic 
matter at its core. The presence of the exotic matter allows for the existence of the non-trivial topology. As exotic 
matter, a massless ghost scalar flcld has been chosen (see the Lagrangian ([T])). The perfect fluid has been modeled 
by a polytropic equation of state For this system, by means of numerical calculations, we have obtained regular 
solutions, which we demonstrated for the parameters 7 = 5/3 and 7 = 1.5 (see Figs. [T] and [21 respectively). 

Our general procedure to obtain the solutions is summarized as follows: we start from the core of the configuration 
(from the throat of the wormhole) with some initial core density of the fluid corresponding to the density at the 
throat. This core density can be chosen either as the density of an ordinary non-relativistic star or, for instance, as 
the density of a neutron star. The condition of regularity of the solutions at the core then requires an appropriate 
choice of the boundary conditions (see Eqs. ([Tl]) -([T5 |) ). The subsequent behavior of the solutions is determined by 
the model's parameters: the polytropic index n and the parameter a which denotes the ratio of the core pressure Pc 
and the core density pc of the fluid. As the dimensionless radial coordinate ^ increases, the density of the matter 9 
decreases, reaching zero at some boundary value = £,b- This value of the radial coordinate can be considered as 
the boundary of the ordinary matter. From the point of view of a distant observer, this is a visible boundary of the 
star. Note, that in this model the scalar field (f> and its derivative 4>' are not yet equal to zero on the boundary of the 
fluid. It requires matching of the internal fiuid solutions with external solutions having a nonzero scalar field energy 
density, that vanishes asymptotically. 

The essential point is that this model differs from ordinary stars by the presence of a tunnel at the star's core. It 
is obvious that, in principle, the matter can move freely through the tunnel. To investigate such motion, one can 
use the approach suggested in Ref. [18[ to describe infinitesimal radial oscillations of the gas within the framework of 



general relativity. Applying this method to the "star plus wormhole" system, we have derived the system of equations 
and describing the infinitesimal radial oscillations of the configuration. For a polytropic equation of state, 
the investigation of this system of equations is presented in Appendix [Bj As in the case of relativistic stars without 
wormhole discussed in Ref. [isj, solving the equations needs a special approach because of the presence of singular 
points in the equations (|B.11[) and (|B.12p . An example of a solution is given in Appendix |B] (see Fig. ^ . 
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Another important question arising when considering models of eompaet objects is the stabihty of such config- 
urations. The question of stabihty of wormholes with massless ghost fields or phantom matter was, for instance, 
investigated in [13, 120l - l23j (see also references therein) . As discussed in [H, [131 j the main difficulty of such studies 
is related to the behaviour of the perturbations near the throat. For the "star plus wormhole" system in Appendix 
[Cl a qualitative analysis of the possibility of the existence of stable solutions outside the fluid was carried out. This 
analysis has shown that the use of a massless ghost scalar field as a source of nontrivial topology allows for the 
existence of an infinite number of oscillating perturbed modes whose energy, nevertheless, remains finite. To reach a 
definite conclusion concerning the stability of the system under consideration with respect to radial perturbations of 
the scalar field, however, the full system must be analyzed. 

The investigations presented in the paper can be considered as preliminary ones. They show the principle possibility 
of the existence of compact astrophysical objects of the type "star plus wormhole". Further investigations can be 
directed towards: (i) a more detailed analysis of such type of configurations with different parameters of the polytropic 
equation of state; (ii) the use of a more realistic equations of state, for example, for the description of the structure 
of neutron stars when it is important to take into account relativistic effects; (iii) the estimation of the energy 
associated with the presence of the oscillations of stellar objects, and possible applications of such estimates for a 
description of various astrophysical events in the Universe (for example, the gamma-ray bursts); (iv) the performance 
of a comparative analysis of the obtained models with the known models of ordinary stars with trivial topology for 
the purpose of identifying the differences, which can lead to a number of observational consequences. 
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Appendix A: Transformation to a new variable 

To perform an analysis of the system "star plus wormhole" which is described by the Lagrangian ([1} and a polytropic 
equation of state (jH), it is more convenient to introduce the new variable 9 related to the density p at a given point 
and the core density pc by 



where n is a constant related to 7, Eq. (|4]), as 71 = 1/(7 — 1). Such variable allows to perform a comparative analysis 
of the model both for relativistic and non-relativistic polytropic stars. Then, according to Ref. nj, we choose the 
static metric in the form 



sm 



(A.l) 



where the metric functions v, A depend only on the radial coordinate r. The corresponding gravitational equations 
will then be 



G\ 



1 








1 ^ 











+ ^ = SttGT/ = SttG 



1 



-A J./2 



(A.2) 
(A.3) 



where the energy-momentum tensor from Eq. ([2|) has been used for obtaining the components Tq and T^. The 
expression for 0'^ can be obtained by integrating the equation on the massless scalar field 



1 d_ 

l—g dx^ 



-99 



= 



in the following form 



n2 



(A.4) 
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where D is an integration constant. Then the system of equations (|A.2[) - (IA.3I) takes the form: 



1 

l2 



One more equation foUows from the conservation law T^.j^ 



- SttG 

- SttG 
0: 



1Z?2 

+ 2—' 



dp 
dr 



1 , dv 



Following Ref. [17| . we now define a new function u by 



u{r) 



2GM 



(1 



2GMu 



(A.5) 
(A.6) 

(A.7) 
(A.8) 



where M is the mass of the configuration within the range < t" < rf,, where rf, is the boundary of the fluid where 
6 = 0. Using this function, Eq. (jA.5[) becomes 



M 



du 
dr 



P 



2 r* ' 



(A.9) 



One can see from this equation that Mu(r) = A/(r) can be interpreted as the total mass of the configuration, 
including contributions both from the fluid and the scalar field, within a sphere of coordinate radius r. For the 
spherically symmetric case without a wormhole (what, in our case, corresponds to absence of the scalar field), one has 
to put it(0) = to avoid a singularity for the mass at the origin [13] ■ This corresponds to the fact that the mass at the 
origin is equal to zero. In the model with the wormhole, there exists a minimal value of r = rg, and correspondingly 
there exists some minimal mass. Therefore we have u{ro) ^ (see Eq. pT|) ). 
Next, using we have from (|A.7p 



with 



2a{n + l)de + (1 + (je)dv = 0, 



(A.IO) 



Pc 



(A.ll) 



where Pc is the pressure at the core of the configuration. Eq. (|A.10[) may be integrated to obtain in terms of ( 



1 



2(n+l) 



(A.12) 



where e"^ is the value of e'^ at the core of the configuration where 9 = 1. The integration constant j^c, corresponding 
to the value of v at the throat, is determined by requiring e'^ to be equal to unity at infinity that corresponds to an 
asymptotically flat space-time. 

Next, let us rewrite Eq. (|A.6p in the following form: by using (|A.9p one can exclude from (|A.6[) the term 



IE: 

2 73 



1 i,,du 
rM — , 



47rr^ dr ' 



and express from (jA.lOP the derivative 



2(t(?i+ 1) de 



l + aO dr' 

Substituting the two last expressions into Eq. (jA.6p . and using the definition for e^^ from (jA.Sp . we find 



- 1- 



2GMu 



1 



1 2cr(n 
r 



I) de 



1 + ae dr 



1 



8ttG 



(A.13) 
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Now we make a change of variables which puts Eqs. (jA.9|) and (jA.13[) in dimensionless form, defining 

1/2 



^ = Ar, u(^) = u(r), where A 



(A.14) 



with A having the dimension of inverse length. (Here we use the dimensionless variable according to the notations 
of Ref. ■ In section IIVI we have already used the letter v for the designation of the three- velocity. But this should 
not lead to confusion.) Using these new variables v{^) and ^, Eqs. (jA.9|) and (jA.13[) take the form 



I -2a{n + l)v/^ de 



1 



a6 di 
1L>2 



r (1 - ad) 



1 dv 



2 e 



1 



-2(ri+l)' 



(A.15) 
(A.16) 



with a new dimensionless integration constant 



D = 



AttGD 

{n+l)a' 



'Pc 



In the absence of the scalar field, i.e. when D = 0, this system reduces to the known system of equations for the 
relativistic polytropic gas of Ref. [131 ■ the non-relativistic limit, i.e. for cr (corresponding to p <C p), the system 
of equations (|A.15I) - (|A.16|) without the scalar field reduces to the Lane-Emden equation 



,d£ 



Appendix B: Oscillations of the system with a polytropic equation of state 



As note above, oscillations for various equations of state of matter creating compact configurations have been con- 
sidered extensively earlier. In particular, for the polytropic equation of state © used in Ref. |18| for the investigation 
of the dynamical instability of a gas sphere, the corresponding perturbed component of the pressure Pp is given by 



hence 



Taking into account that 



we have 



O'p + nfe. 



Pp = uPcOq ^6p , 



Pp+Pp = PcOq 



Next, 



Po+Po^ PcOo{l+<jeo) 



- + a{n + l) 



p'o+p'o = pM 



+ cr(n+ 1) 



Using the obtained expressions, equations (|38p - pn|) can be rewritten as 



,-Ao J,' 



Xp ^ -SnGe^^r [p.e'Sil + aOoK 
npctiQ z r 



-Pc(^a 



,-An 



A,] 



p r 



(B.l) 
(B.2) 
(B.3) 
(B.4) 



(B.5) 

(B.6) 
(B.7) 
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FIG. 3: The dimensionless displacement tp obtained by solving the system of equations HB.11|) and (|B.12|) for the case n = 1.5. 
For these curves, the values of oj^ are the following: for the top dashed line cj+ = 0.177, for the bottom dashed line lu'L = 0.175. 
The solid line corresponds to oj^ — 0.1753 and A ~ 0.001. The boundary of the fluid is located at ~ 1.79. 



These equations can be rewritten through the dimensionless variables from ()A.14[) and new dimensionless variables (p 
for the scalar field and tjj for the displacement 



/Pc 



V, c 



A' 



Taking into account that (see Eqs. (|A.8p and (|A.14p ) 



Vq 



l-2a{n+l)j 



^, ^ 2a(n + l)K-i.o/g) 
" ^-2(7{n + l)vo ' 



and the expression for i^q from (|A.12[) 



1 + 0- 



1 2(71+1) 



1 + (t6 



1 + (t6 



we find from ((R5)) - (|BJ)) : 

Ap = -2a(n + 1)6^"^ K(l + ^^o)V^ + e-^Vo<Pp] , 

+ e-^" [Kv'o-V>o] ^p}, 



(B.8) 

(B.9) 
(B.IO) 



where "prime" denotes now differentiation with respect to the dimensionless variable ^. 

Using the obtained dimensionless expressions, Eq. (|42p takes the form (here a new dimensionless frequency Co is 
introduced, u) ~ '-o j A\ we further drop the bar for simplicity) 



91 + 



+ cr(n+ 1) 



(B.ll) 
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The perturbation equation for the scalar field (|43)) is then given by 

+ \ - ^'o + ^^e^»-^Vp = 0. (B.12) 

So we have two linear second-order ordinary differential equations (jB.lip and (jB.12[) for the dimensionless pertur- 
bations of the scalar field (fp and the displacement ip. 

As an example, let us demonstrate the solution of this system for the case of n = 1.5 whose static solutions have 
been found in section lllll (see Fig. [1]). The numerical calculations show that for this case the solution of the system 
of equations (jB.lip and (|B.12p has a singularity at the boundary of the fluid when 6* — > 0. For some value of the 
frequency cj+ the solution for the displacement ip tends to -l-oo, and for some value u>- to — oo (see Fig. [3|). One 
expects that there exists some value of uj in the range between and w_ at which the solution becomes regular. In 
order to find such solutions, it is possible to use one of the methods described in Ref. [l^. One such method consists 
in the following procedure: determining the eigenfunctions accurately everywhere, one can integrate simultaneously 
from the throat and from the boundary of the fluid, and require matching of the logarithmic derivative of ip at an 
interior point. On the other hand, it is possible to use the shooting method when one can try to find (using a method 
of step-by-step approximation) an eigenvalue of w at which the eigenfunctions are regular. Using such an approach, 
and choosing w, one can try to achieve that the solution will be approaching closer and closer to the boundary of the 
fiuid until it will not reach some critical value £,crit at which the solution diverges. The width A = (^j, — S^crit) is 
determined by: (i) the structure of equations (jB.lip and (jB.12p : (ii) the accuracy of the numerical method used; (iii) 
the machine precision. Using the NDSolve routine from Mathematical we have reduced A to 0.001 at ux^ ~ 0.1753. 
The solution for this value of is shown in Fig. [3] by the solid line. 



Appendix C: Stability analysis 



In this section we study the dynamical stability of the above static solutions under linear perturbations. Due to 
considerable technical challenges in performing the stability analysis for the system as a whole, we just present a linear 
stability analysis of the external solution outside the fiuid. Restricting ourselves to this region, makes the calculations 
much simpler, while it is still possible to demonstrate some general properties concerning the stability of the system 
under consideration. 

Outside the fluid, the static solutions are given by the system of equations (|^ ^ - (|^^ . The corresponding time- 
dependent equations outside the fluid can be obtained from equations (j^ - (j?T|) by putting p = p = Q. The time- 
dependent equation for the scalar field can be written as follows: 



-A 



J," 

9 



= 0. 



(C.l) 



We perturb the solutions of this system by expanding the metric functions and the scalar field function to first order 
as follows: 



A = Ao(r) + Ai(r) cos wt, v = i>Q[r) + vi{r) cos ujt, (p ^ (f)o[r) + (f)i{r) 



cos Lot 



(C.2) 



The index indicates the static background solutions, and the index 1 refers to perturbations. Substituting these 
expressions into the Einstein equations (j26p - (j28p . one can find the following expressions for the metric perturbations: 



Ai = -87rG0o<? 



and 



.-a;-^[2- 



r{v'^-K)] cPi 



(C.3) 



(C.4) 



Then, substituting the perturbations (jC.2p and expressions (jC.3P and (jC.4p into Eq. (jC.ip . one can obtain the 
following equation for (pi rewritten in the dimensionless variables used in the previous sections: 



with the potential given by 



VoiO - ^^^^ + + Wo [2 + K)] > 



(C.5) 



(C.6) 
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u.uu i- I 

rij 5 10 15 20 r| 

FIG. 4: The form of the potential (|C.9|) for the case of n = 1.5, a — 0.2. The potential is always positive for these parameters. 
The left boundary is situated at the edge of the fluid at rjt = £,b ~ 1.79. Asymptotically, aX rj ^ oo, the potential goes to zero. 



and io = Lj/A. Introducing the new independent variable rj 



= cxp 



1 



(Ao - vq) 



one can rewrite Eq. (jC.5[) in a Schrodinger-likc form 



where 



(C.7) 



(C.8) 



(C.9) 



Now the consideration of the stability of the system reduces to a study of Eq. (|C.8P . Then if the energy- like eigenvalues 
are positive, the solution is stable w.r.t. the perturbations addressed, due to the regularity of the perturbations (see 



the expressions for the perturbations (|C.2[) ). Another important point, providing regularity, is the localization of the 
eigenfunction (pi in a bounded region of space. The qualitative behavior of the solutions of Eq. (jC.SP can be estimated 
from the form of the potential (jC.9p . As an example, in Fig. |4]the form of this potential for the case considered in 
Appendix |B] with n = 1.5, a = 0.2 (see Fig. [3]) is presented. One can see from Fig.|3]that, using the above-mentioned 
parameters, the potential (|C.9|) looks like a wall near the boundary of the fluid at ^ = ^f, (the variable r] is chosen in 
such a way as to provide rjb ~ S,b at the boundary of the fluid) . In quantum mechanical language, this will correspond 
to motion of a particle in the field of the potential (jC.9|) with an oscillating eigenfunction (fi. The absence of a well 
in the potential (|C.9p inhibits the existence of a discrete spectrum of energy levels of the system (|C.8P - (jC.9|) . 
Let us now estimate energy associated with these oscillations. The energy density of the scalar field is 



Then using the expressions for perturbations from (|C.2|) . one can obtain the perturbed part of (jClOp 
form 



(C.IO) 



in the 









[;( 







cos wt, 



or in dimensionless variables 



cos WT, 



(C.ll) 



where t = At is the dimensionless time variable, and a new dimensionless constant 
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is introduced. Then the proper energy of perturbations is defined by the integral 

E = 

or in dimensionless variables 



E — An cos UJT 



(C.12) 



where E = E / E^. \s the dimensionless energy measured in units of some characteristic energy Ec = pc/A^- To show 
that the oscillations have finite energy, we can perform an analytical estimation of the asymptotic value of this integral. 
To do this, let us find asymptotic expressions for the functions <y9g, Aq, Ai and (jyi appearing in the integral. It follows 
from equation (|22p that asymptotically 



di 



152 



52 
2^ 



(C.13) 



Then from p7)) we have 



Next, equation gives 



1 + a{n + 1) 



HI 

e 



The asymptotic expression for 0i can be found as follows: one can see from the form of the potential (|C.9|) that it 
approaches zero asymptotically. Correspondingly, equation (|C.8p describes motion of a free particle, and its solution 
can be presented in the form 

(f)i « a cos (oi^ — /3), 

where a, /3 are some arbitrary constants. Then the expression for Ai from (|C.3[) takes the form 

^aeos (lj^ — /?) 



Ai « -8nGy/p;D- 



e 



Using all these asymptotic expressions, we have the following asymptotic form of the integral from (jC.12p : 

sin (w^ — /3) 



E « AiraDDu cos 



-dL 



(C.14) 



where the initial integration point is chosen at that place where the potential (|C.9[) is close to zero. Since this 
integral is converging, the proper energy of the oscillations remains finite. Thus, taking into account that the value 
is always positive, one can conclude that the external solution is stable against these linear perturbations. 
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